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000000 TTTTTTITTT SSSSSSSS PPPPPPPP 000000 WW WW HH HH 
000000 TITTITTTTT SSSSSSSS_ PPPPPPPP 000000 bald Ww HH HH 

00 00 TT SS PP PP 00 00 WW WW HH HH 

00 00 TT SS PP PP 00 00 WW WW HH HH 

00 00 TT S$ PP PP 00 00 WW Ww HH HH 

00 00 TT SS PP PP 00 00 WW WW HH HH 

00 00 TT SSSSSS PPPPPPPP 00 0O WW WwW HHHHHHHHHH 

00 00 TT SSSSSS PPPPPPPP 00 00 WW WW HHHHHHHHHH 

00 00 TT SS PP 00 00 WW WW WwW HH HH JJ 

00 00 TT SS PP 00 OO WW WW Ww HH HH JJ 

00 00 TT SS PP 00 00 WWWW WWWWw HH HH JJ 

00 00 TT ss PP 00 00 WWhW WWWw HH HH JJ 

TT SSSSSSSS PP 000000 tel wl Ww HH HH JJJJIJ 

000000 TT SSSSSSSS PP 000000 tlw Ww HH HH JJJJIJ 

Lt. III SSSSSSSS 

LL HI] SSSSSSSS 

LL I] SS 

LL I] SS 

LL I] SS 

LL I] SS 

LL I] SSSSSS 

LL I] SSSSSS 

LL 1] SS 

LL I] SS 

LL I] SS 

LL I] SS 
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9) et DEEL AR ; Detailed Current Edit History 
(4) 107 OTs$ PCO eR - REAL*16 ** INTEGER*®4 
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= REAL*16 ** INTEGER®4 power ro nine 16-SEP-1984 AX/VMS Ma v04-00 P 1 
— ’ -SEP-1984 2 94 9B: 33 ERTHRTL. SREJOTSPOWHS MAR; 1 = (1) 


OWHJ = REAL*16 ** INTEGER*®4 Royer routine 
4/ ; File: OTSPOWHJ.MAR Edit: SBL1004 
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—MARABARARAARAARARAAALAAAAALALALE AAAS RESE SLES ESSERE RRR ERE EERE ESRC E EERE RRR R SESS 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED 


** ® 
*® ® 
:* ® 
‘® ® 
:* ft 
i THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
i ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
te INCLUSION OF THE ABOVE COPYRIGHT NOTICE, THIS SOFTWARE OR ANY OTHER © 
;* COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
;* OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
s*# TRANSFERRED. * 
** ® 
it THE INFORMATION IN THIS SOFTWARE 1S SUBJECT TO CHANGE WITHOUT NOTICE 
;* AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
it CORPORATION. . 
:* ® 
** *® 
*® ® 
:* ® 
:* * 
*® ® 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


(ARRRRASARASAAALEARSASELESEL ESAS EASELS SEES ERAS RRR RSE S SEER ERR R ESE S ESR R EERE 


s FACILITY: Language support Library - user callable 
; ABSTRACT: 


REAL*16 base to INTEGER*4 power. 

Floating overflow and underflow can occur. 

Undefined exponentation can occur if base is 0 and power 
is 0 or negative. 


VERSION: 1 


HISTORY: 
; AUTHOR: 
Steven B. Lionel, 11-Jun-79: Version 1 
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-SBTTL HISTORY i Detailed Current Edit History 
4; pb $ History for Version 1 of OTSSPOWHJ 
0 5 ; 1-001 = Adapted from OTS$POWDJ version 1-001 SBL 01-Jun-79 
4 $ 3 I~856 - Recast some comments. JBS 30- 979 
00 ; 1-005 = Add handlers to catch SS$_FLTOVF and SS$_FLTDIV, _ sf nal 
000 8 : MTHS_FLOOVEMAT or MTHS_FLOUNDMAT instead, dependin he context. 
00 59 : Also disable Iv ang sonan e BLBS/ASHL at 5x6 R to BBSC/ROTL for 
0 60 ; uniformity Sith OTS W 26-Feb-1 
444 o3 ; 1-004 = Use word-relative a ei pol handler address. SBL 9-Sep-1981 


sap temeanreaithcianiaasiteatstnit — —+ 


3 
OTS$POWHJ - REAL*16 ** INTEGER®4 power routine 16=SEP=1984 02:01:09 VAX/VMS Macro V04-00 Page OTS’ 
a 806 DECLARATIONS 2 mets fi 9950858 MTHRTL.SRCJOTSPOWHJ .MAR; 1 ° 3) 1-8 
rt: -SBTTL DECLARATIONS 
96 : 
67 ; INCLUDE FILES: 
OB 
3 ; 
4 3; EXTERNAL SYMBOLS: 
000 * 
0 74 -DSABL GBL 
8 75 sEXTRN MTHSK_UNDEXP, MTHSK_FLOOVEMAT, MTHSK_FLOUNDMAT 
4 ~-EXTRN MTHS$SSIG ; Math error routine 
909 A .EXTRN SSS$_FLTOVF, SSS_FLTOVF_F, SSS$_FLTDIV, SS$_FLTDIV_F, SS$_CONTINUE 
0000 79 ; MACROS 
0000 80 ; 
0000 81 SCHF DEF ; Define condition handler symbols. 
0000 Hf SSFDEF 3; Define stack frame symbols. 
0000 SPSLDEF ; Define program status longword 
0000 He 3 symbols. 
0000 5 
0000 86 ; 
0000 87 ; EQUATED SYMBOLS: 
0000 88 ; 
00000004 0000 89 base = 4 3; base input formal - by-value 
00000014 0000 90 exp = 20 ; exponent input formal - by-value 
0000 91 ; Note: quad loa tne by-value violates 
0000 4 ; calling standard, but ok since this 
0000 & ; routine is a code support routine (OTS$) 
0000 94 
0000 95 ; 
44 96 ; OWN STORAGE: 
000 97 ; 
0000 98 
0000 99 ; 
000 100 ; PSECT DECLARATIONS: 
000 101; 
000 13 
00000000 10 -PSECT _OTSSCODE PIC,SHR,LONG,EXE,NOWRT 
444 H+ ; program section for OTS$ code 


r 4 
POWHJ_R3 = REAL*16 ** INTEGER®4 6-SEP-19 MTHRTL.SRCJOTSPOWHJ .MAR; 1 (4) 


} 
3 | 

OTSSPOWHJ = REAL*16 ** INTEGER*4 power vevites 16-SEP-1984 :01:09 VAX/VMS Macro V04-00 P 4 
80. ores 4 gY 949858 ( — 
~SBTTL OTSSPOWHJ_R3 = REAL*16 ** INTEGER*4 | 


;++ 
FUNCTIONAL DESCRIPTION: 


REAL*16 result = REAL*16 base ** signed longword exponent 
The REAL*16 result is given by: 


base exponent result 


any >0 product (base * 2**i) where i is each 
non-zero bit position in exponent 


(Pe 


v AUv 
o ooo 


product (base * 2**j) 
is each non-zero bit position 
ponent: 


—- 


ned exponentation 


Al 
oo 
AA 


0 
0 
e 
d 
0 / product (base * 2**j) 

ere i is each non-zero bit position 
in ‘exponent: 


3 Floating overflow can occur on either of the two MULH's. If this 

; aapoens when the exponent is less than zero, the except ren is caught by 
3 a local condition handler named EXC_HNDLR_UNDER, which sets the result 
3 to 0.0 and either signals MTH$_FLOUNDMAT Tif FU is enabled in the 

3 caller's PSW) or continues at POWHJX. If it happens when the exponent 
3 is greater than zero, the exception is cong by a local condition 

3 handler named EXC_HNDLR_OVER, which sets the result to the reserved 

: operand (-0.0) and signals MTHS_FLOOVEMAT. 


Floating overflow and floating divide by zero can occur on the DIVH. 
These exceptions are caught by EXC_HNDLR_ OVER, which sets the result to 
the reserved operand (-0.0) and signals ATHS_FLOOVEMAT. 


Undefined exponentiation occurs if base is 0 and 
exponent is 0 or negative. 


CALLING SEQUENCE: 
Power.wh.v = OTSSPOWHJ_R3 (base.rh.v, exponent.rl.v) 

INPUT PARAMETERS: 

; base - 
exponent - 


IMPLICIT INPUTS: 
The setting of FU in the caller's PSW. 


OUTPUT PARAMETERS: 
NONE 


6 base 
R 


} 
= 0 1.0 | 
. 8 Undefined exponentation 
*4 exponent 


SOOOCOOCOOOOCOSOOSOSOSOOOOOSOO OOOO OOOOSCOOOCOOOOOOOO OOOO OOOOoOO 
SOOCOCOCOCOSOSOOSOSOSSOSSOSOSOSOOCOSOSCOCOSCSOOOSOOSOOOSOOOOOOOOSOOOOOOOOOOOOO 
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DPDPRDPRD PP Be BS BBE BE ENN WINN WIAIDIPINIPOPONNINID 2 2 OOOO OOOO 
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; IMPLICIT OUTPUTS: 
NONE 


fe 
2 


OTSSPOWHJ - R 
1 OTS 


v> 
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) 16-SEP-1984 96:01:08 palyas 


r routin v04-00 P 
NTEGER 4 6-SEP=1984 11:28:28 CMTHRT S$ sat 


acro 
RCJOTSPOWHJ .MAR; 1 


a 
Wwe 


FUNCTION VALUE: 
WARNING !! 
REAL*16 result in registers RO through R3. 
This is a serious violation of the VAX calling standard and can 
kit programs if misused. Beware! 
SIDE EFFECTS: 
Destroys registers R2-R3! 
Signals MTHS _FLOOVEMAT if ge te overflow occurs on either of the two 
_ $s & wien, on rae > 0, f floating overflow or divide by zero 
Signals & nts F LOUNDMAT Hy Py goby over flow eo on either of the two 


ULH's en ex xponent aller has FU enable 
SIGNALS MTHS UNDE = ‘, UNDER INED EXPONENTATION® 5 if 


base is 0 and re is 0 or negative. 


-ENTRY OTSS$POWHJ_R3, “M<R4, RS, RO, R7, RB> 
on't save regtacers R2 and R3! 


3 Disable sneeorr overflow. (Occurs on 
; maximum negative exponent.) 
6D OO6F'CF 9E MOVAB W*EXC_HNDLR_OVER, (FP) ; Translate exceptions to 
3; MTH$_FLOOVEMAT. 
MOVH #1, RO + RO-R3 = initial result 
MOVH base(AP), R4 3 ReaR? = base 
MOVL exp(AP), R8 ; RB = exponent 
BGTR EXPGTR ; pranch if exponent > 0 
MOVAB B*EXC_HNDLR_UNDER, (FP) ; Translate exceptions to 
3 mts _FLOUNDMAT. 
TSTH R4 3; test base 
BEQL UNDEF INED 3 undefined O**0 or 0**(-n) 
MNEGL R8, RB : = fexponenti 


BEQL POWHJX if exponent is 0, return RO = 1.0 


a¢ 
; Exponent is > 0 or (exponent is =< 0 and base is not = 0 -- use ‘iexponent:) 


“WH DOODOOOCOOCO OOO OOOO OOH 0000000000 00 09 09 09 SI NINN SINS NS NNN AOA Oo 
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SOOWONO UNS Vt ~* . DPBNA MAE WN OS OD NA UE WIN @ 0 OO NAMES WN 0 OONOU ES WN OOONOUNS 


and signal SS$_FLTOVF or S s FLT OVF _F. 
branch if next bit in ‘exponént: is” 


"XPGTR: BBSC #0, R8 wage, ; branch if ‘exponent i is odd 
SQUAR: ROTL al, ar’: ; RB =_‘exponent:/2 
SQUAR1: MULH2 3; R4-R7 = current power of beee 
; floating overflow will trap or fault 


BLBC R8, SQUAR 


; Here when bit i of ‘exponent: is a } 


> Partial result = partial result * base * 2**j) 


AAAAA. A. AARINININININININ 4 9 9 PF FH QOOCOOCOOCOOCOCOOCOSOOSOOCOOOOOSOOOOOOOOSOOOSOOOoO = = 
BPE MO DHL LL LL HSMN PS POLO FNVNVININVOOCOCOCOCCOOCOCCOOCOCOOCOOCOOOOCOO wae 
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OTSSPOWHJ - REAL*16 ** INTEGER®4 power routine 16-SEP-1984 01:09 VAX/VMS Macro V04-00 Page 6 
9 80. OTSSPOWHI_RS - REAL*16 + INTEGER*4 ats y et BE: EMTHRIL. SREIOTSPOWHY MAR: 1 . (4) 
PARTIAL: 
50 54 64FD MULH2 R4, RO : arte = new partial result 
58 58 FF BF 13 ASHL #1, RB, RB : RB = lexponent! 
cE 1 BNEQ SQUAR1 ; loopback if more exponent bits are 1 
14 AC D5 TSTL exp (AP) 3; test sign of exponent | 
9 14 BGTR POWHJX : if exponent > 0, return RO 
6D O6F'AF  9E MOVAB B*EXC_HNDLR_OVER, (FP) ; Translate exceptions to | 
3; MTH$_FLOOVEMAT. 
50 08 £50 67FD DIVHS RO, #1, RO 3 RO-R3 = 1.0/result 
04 POWHJX: RET 3 return, result in RO-R3 


° 


UNDEF INED: 
> Ri FF Ff ASHQ 


52 
OO'8F 9A 


CLRQ 
7E MOVZBL 
00000000'GF 01 FB CALLS 


MIMI WANN 


OOCCCCCCOCO BB NOL MMMMMmMmo Cocos rye MOC Se & 


RET 


34 
; The following 
3 when exp < 0) 


PAAASOCAOC AOC w SPP EES 


3; underflow of t 
; MTH$_FLOUNDMAT. 


3¢ 
3 Undefined exponentation error - 0**0 or 0**(-n) 


#15, #1, RO 


R 
#MTHSK_UNDEXP, =(SP) 
#1, G*ATHS$SIGNAL 


RO-R3 = reserved floating operand 

second quadword 

Indicate undefined exponentiation. 

convert to 32-bit condition code 

and SIGNAL MTH$_UNDEXP 

Note: 2nd arg not needed since no JSB OTS$ 
is possible. 

return 


handler is established to process exceptions which imply 
he final result (floating overflow in either of the two MULH's 
On the occurrence of such an exception, the handler signals 


OOSCOCCCOCOOCOOOCOOOOCOOSOOCSOOOOOCOOOOO OOOO OOOOOOOOOOOOO 
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EXC_HNDLR_UNDER: 
-WORD “M<R4, RS, RO> 3 Entry mask. 
§ ; Don't save registers R2 and R3, 
28 «10 BSBB SETUP ; Set up RO:R5 and identify condition. 
64 3; Return only if FLTOVF or FLTDIV. 
1E 04 AG «06 «ET =—(0064 BBC #PSLSV_FU, SFS$W_SAVE_PSW(R4), CON_U 
69 ; Branch Tf caller has not enabled FU. 
56 00'8F 9A 69 MOVZBL #MTHSK_FLOUNDMAT, R6 ; Report MTH$_FLOUNDMAT, not SSS$_FLTOVF. 
ao er BRB DO_SIG 
F 3+ 
oF ; The following handler is established to process exceptions which impl | 
6F ; overflow of the final result (floating overflow in either of the two MULH's 
O6F 3; when exo > 0, floating overflow in the DIVH, or phase ine divide by zero in the 
O6F : DIVH). On the occurrence of such an exception, the handler signals 
O6F ; MTH$_FLOOVEMAT. 
O6F i- 
O6F 
6F EXC_HNDLR_ OVER: 
0070 of -BORD “M<R4, RS, RO> ; Entry mask 
1 ; Don't save registers Re and R3, 
1¢ 6.10 «=90071 BSBB SETUP 3; Set up RO:R5 and identify condition, 
073 3; Return only if FLTOVF or FLIDIV. 


4 
TSSPOWHJ AL*16 ** INTEGER®4 power noun lon y is AX/VMS Macro V04-00 Page 
I 856 OTSEPOUHI_AS ORES eTG Oo ENTE ERC “1386 7:5b:58 EMTaNte caedoreocucy.mar;1 29 
50 01 OU 7 7 #15, ; Make the default result -0.0. 
56 00'S 9A rt hi mOVIBL PHThisK ELOOVEMAT, R6 ; Report aTHs. FLOOVEMAT, not $$$_ FLTxxx. 
10 ASOD 78 1 DO_SIG: PUSHL SFSL_SAVE_PC(R4) ; Report caller's PC, not exception PC. 
5 DD 7 g PUSHL R : Report MTH$ xxx, not SS$_xxx. 
00000000 ' GF ? FB CALLS » G*MTHSS$SIGNAL ; Signal fre Condition. 
OC AS 7D 4 CON.U: MmOVa RO, CHFSL_MCH_SAVRO(RS) ; if wy gued jfestece RO and R1. 
8 5 ; and R3 will retain their values. 
50 00° pO 0088 § MOVL S“#SS$_CONTINUE, RO ; eet leen wilt the original exception. 
04 Ht a DO_RET: RET ; Exit from handler. 
F 3+ 
38 ; 83 ; Common etye routine for handlers. Returns normally if exception was FLTOVF, 
O8F 91 ; FLTOVF_F Tolv. or FLTDIV_F. If the exception was nangthing else, it 
st 38 ; executes a RET feusing 8 an exit from the handler with which is 
O8F 93 ; equivalent to §s$ ngs} NAL. In the gees of a normal return (FLTOVF, FLTOVF_F, 
OO8F 94 ; FLTDIV, or FLTDIV_F) it aie up RO:R5 as follows: 
OO08F 95 ; RO/R3: O 
OO8F 96 ; : address of establisher's frame 
OO8F 97 ; RS: address of mechanism array 
Bh Bap 
50 7C¢ SOB 00 SETUP: CLRQ RO ; Clear first half of default result. 
54 04 AC 7D 0091 01 mova CHFSL_SIGARGLST(AP), R4 ; R4 = address of signal array 
0095 ¢ R5 = , oddress of mechanism array 
0000 * 8F 04 AS «BI Ssé0095 0 CMPW CHFSL_SIG_NAME(R4), #SSS_ gj bge 
0098 04 ; Was it a floating overflow trap? 
18 13 0098 05 BEQL DO_RSB Branch if yes. 
0000* BF 04 AG «Bl OO 06 CMPW CHFSL_SIG_NAME(R4), #SSS_ gyn F 
OOA 07 ; Or a Floating overflow fault? 
10 13 OOA 08 BEQL DO_RSB Branch if yes. 
0000 ' 8F 04 AS =©B1 OCOQOA 09 CMPW CHFS$L_SIG_NAME(R4), #SS$_ FLIDIV 
QOAB 10 ; Or a floating divide by zero trap? 
08 13 OOAB 11 BEQL DO_RSB Branch if yes. 
0000 * 8F 04 AS «Bl SOCOQOA i CMPW CHFSL _SIG_NAME(R4), #SSS$_ FLIDIV F 
008 1 ; Or a Floating divide by zero fault? 
09 =«12 set 3 BNEQ DO_RET : mene of the above: return from handler 
08 AS) 9S AF 9— 0085 16 DO_RSB: MOVAB B*POWHJX, CHFSL_SIG _NAME +4 (R4) 
OOBA 17 ; Change return PC to POWHJX. 
54 04 AS 0 OOBA 18 MOVL CHFSL_MCH_FRAME(R5), R4 ; R4 = address of establisher's frame 
Beek 19 CLRQ R2 : Clear second half of default result. 
Beet 31 RSB ; Return. 
b0e1 2 END 


(4) 
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OTSSPOWHJ - REAL®16 ** INTEGER®4 power ro ste i 4 AX/VMS -00 P TS$ 
Symbol table pew ire sett 4: 33; 3 YatANth L.S E51 ebouny MAR: 1 _ d) abl 
BAS 


ron 
Qzr 


2Bz8 


mmmm 


| 


P~ 
PGTR 2 
MTH$$S1GNAL eeeerere 
MIHSK_FLOOVEMAT eeeeeeee 
MTHSK~ UNDE XP eeerener 
OTSSPOWHJ _R3 900 RG | 
PARTIAL 00034 R 
Jx o'9 0 & 
PSL$V_FU = 
TUP OB 02 
SFSL_ = 10 
SFSW SAVE PSW = 
SQUA 00000028 R f3 
SQUAR1 000000. > R 0 
$S$$_CONT INUE eeererne x 00 
SSS"FLIDIV terereee xX 00 | 
SSS$"FLIDIV_F eereeeee xX 00 | 
SS$_FLTOVF eeeeeere x 00 | 
SS$-FLTOVF_F teverers§ xX 00 
UNDEF INED OOO0004E R 02 
$ewecceccesoocccn + 
: :_Psect synopsis ! 
eee eee ee ee ee + 
PSECT name Allocation PSECT No. Attributes 
» 83s 00000000 ( 8. 0.) NOPIC USR CON ABS LCL oo NOEXE NORD ag NOVEC BYTE | 
SABSS 00000000 < 0.) Of 3°} NOPIC USR CON ABS ~ LCL NOSHR EXE RD T NOVEC BYTE 
_OTSSCODE 000000C1 ( 19 eh. wen PIC USR CON REL LCL SHR EXE RD NOURT NOVEC LONG 
' Performance indicators ! 
we wre ew wm wees eK TH eewr en ewe em ewe Be | 
Phase Page faults CPU Time Elapsed Time 
- | 


Seabol table output 
Psect synopsis output 
Cross-reference output 


Initialization 3 0:00: 8.9 0:00:01.5 

Command processing 11 : ’ é : if 9 

Syabol table sort 7 00: ; :00: 9 
2 t ; 00:09: 


& | 
GNP Acro ma satiety 7 MATE #7 IMTEGENG over routle’ TEEEIEE 7:81:08 HANATE SESE ABIG mann POOP cs] | 


Assembler run totals 367 00:00:03.65  00:00:17.23 | 


The working set Limit was 1200 peges. 
1 bytes (18 pages) of virtual memory were used to eurter the intermediate gee. 

here were 10 pages of symbol table we al a dt hold 06 non-local ie local symbols. 
source Lines were read in Pass 1, produc object records in Pass 

11 pages of virtual memory were ea” to defines 10 macros. 


teaser ecases eee scam anrnresas > 


! Macro Library statistics ' 


Sunanmenan eee e sean een en oc e + 


Macro Library name Macros defined 


“$255$DUA28: CSYSLIBISTARLET. MLB;2 6 
148 GETS were required to define 6 macros. 


There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWHJ/OBJ=OBJ$:OTSPOWHJ MSRC$:MTHJACKET/UPDATE=(ENHS$:MTHJACKET) #MSRC 


| 0265 AH-BT13A-SE ) NT CORPORATION 
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